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A method for organizing next-to-leading order QCD calculations using a veto which enforces 
the cancellations between virtual and real emission diagrams is applied to hadronic collisions. The 
method employs phase space slicing with the slicing parameter determined dynamically event-by- 
event. It allows for the generation of unweighted events and can be consistently merged with a 
^ . parton shower. The end product is more intuitive for the end user, as it is probabilistic, and can 

^ ' be easily interfaced to general purpose showering and hadronization programs to obtain a complete 

, event description suitable for experimental analyses. As an example an event generator for the 

' process pp^ ^ Z + X at NLO is presented and interfaced consistently to the PYTHIA shower and 

hadronization package. 
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I. INTRODUCTION 

Computer simulations of higher order quantum chromodynamics (QCD) corrections in high energy physics collisions 
typically rely on one of two methods: the parton shower ||l|,|| or numerical integration of next-to- leading order (NLO) 
distributions. Both methods have proven extremely useful, but each has a limited — though complementary — region 
1 '' of applicability. 

Programs which employ the parton shower approach, such as PYTHIA HERWIG Q|, and ISA JET have enjoyed 
^ ■ widespread use by experimentalists. These programs begin with a leading order hard subprocess. Higher order effects 
are added by evolving the event using the parton shower, which provides an all orders description of parton emissions 
valid in the soft and coUinear regions, but is not accurate for well-separated particles. The partons are then grouped 
together into colorless hadrons, resonances are decayed, and the underlying structure of the event (beam remnants, 
multiple interactions, etc.) is added. These programs are able to describe the exclusive structure of the event, and so 
\ provide a useful input for subsequent detector simulation. 

NLO integration programs go one order beyond in the prediction of the cross section, have a reduced dependence 
on arbitrary scale choices, and provide a good description of hard well-separated emissions. They are able to predict 
distributions of events, but are unable to produce events with the frequency predicted by the theory (because the 
cancellations between Feynman diagrams are usually achieved by allowing a fraction of phase space points to have 
negative probability). These programs are commonly used by experimental collaborations for determining fc- factors, 
which are often assumed to be constant, and are employed to correct the normalization of distributions from leading 
order event generators. Lacking individual events to evolve further, it is difficult to add subsequent event features 
like hadronization or the underlying event. This has meant that these programs are excellent theoretical tools 
for predictions of distributions, but limits their usefulness for producing events to be simulated in the detector 
environment. The primary aim of the method presented in this paper is to improve the usefulness of NLO calculations 
for experimental applications by interpreting the result in a manner which is well suited for interface to showering 
and hadronization generators and subsequent detector simulation. 

The pursuit of techniques for combining NLO calculations with the parton shower is a natural direction for the 
evolution of event generators. Two primary challenges stand in the way: overlaps between phase space volumes of 
differing dimensionality need to be accounted for in a manner which does not double-count or neglect any region, 
and the result needs to be interpretable in a probabilistic way (probabilities should be everywhere positive definite). 
Significant advances towards the resolution of the former challenge have been achieved by Collins ||] using a subtractive 
approach, however negative weighted events are still produced. 
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The issue of unweighted event generation has been addressed by the author for the special case of diboson production 
in Ref. Q], then further developed to include a consistent merging of the parton shower in Ref. ||]. In the present 
study, these techniques are generalized and a veto method proposed by Potter |^ is incorporated. This allows for 
the organization of NLO event generators in an elegant and simple manner which competes with leading order event 
generators in terms of efficiency and computer time. 

To illustrate the methods discussed in this paper, an event generator for the process pp^ + X is constructed 

and used to generate event distributions relevant to the Fermilab Tevatron Collider and the CERN Large Hadron 
Collider (LHC). A precise knowledge of the Drell-Yan lepton pair production process is particularly important at 
hadron colliders. It will be used to probe new physics (e.g. large extra dimensions, extra neutral gauge bosons), 
perform precision measurements of electroweak parameters, constrain the parton density functions, and calibrate the 
detector. The latter is perhaps most important to the physics program, because it means our knowledge of this process 
will feed into the systematic errors for most physics measurements: the lepton energy and momentum scale may be 
calibrated in situ with Z'^ — )■ events, the jet energy scale may be determined using events with a Z'^ decaying to 
leptons recoiling against a high transverse momentum jet, and the Drell-Yan event rate may be used to determine the 
absolute luminosity. The NLO QCD corrections to the process have been available for some time Recently the 
complete 0{aQED) corrections have been calculated pl| , and will play an important role for precision measurements. 

In the following section, background information relevant to the NLO calculation is presented. A description of the 
phase space veto method, numerical results, and shower evolution are presented in the next sections. 



II. BACKGROUND 



The NLO cross section receives contributions from the square of the Born graphs, the interference of the Born 
graphs with the one- loop graphs, and the square of the real emission graphs which contain an extra colored parton in 
the final state, 

-^NLO = Born Born one loop + Aureal omission- (1) 

The second and third terms of Eq. |l| diverge when treated separately, and so numerical integration calculations employ 
a regularization scheme which effectively combines pieces of these terms to obtain finite results everywhere in phase 
space. 

Commonly used schemes include the "subtraction method" [|l2|, "dipole method" and "phase space slicing" 
(PSS) methods |14-lJ]. All methods give identical results when used appropriately. For the purposes of this study. 



features of the PSS methods are convenient. Variations of the PSS method include the "two parameter PSS" (see 
Ref. for an accessible review), "one parameter PSS" [|6|, and "Smin-slicing" ]l7| , p^ . 

To illustrate the PSS method, consider pp^ ^ Z'^ + X ai NLO. The Feynman graphs are presented in Figure [|. 
The first term of Eq. |l| is leading order (aqED) and contains n particles in the final state. The phase space volume 
<i>n which defines the configuration of the four-vectors is referred to as n-body and specified by 4 degrees of freedom: 
the Z° mass, Z'^ boost, and two Z^ decay angles. The second term of Eq. || is first order (aqEoas) and is also 
described by n-body kinematics. The third term of Eq. |l| is also first order (q;qedQ!s) and the final state contains 
the vector-boson (or its decay products) and a colored parton (e.g. Z^g, Z'^q, Z^q) and is described by (n-l-l)-body 
kinematics with 7 degrees of freedom: the system mass, system boost, Z'^ mass, two Z^ production angles, and two 
Z^ decay angles. 

For a particular choice of the n-body kinematics, the phase space $+1 which specifies the kinematics of the real 
emission is a plane in u = (p2 — Pj)^ vs. t — (pi — pj)^ space, shown in Figure 0, where pi and p2 are the four- 
momenta of the massless colliding partons, and pj is the massless colored emission (the azimuthal degree of freedom 
is unimportant and not shown). The n-body kinematics occupy a point at the origin of this plane. The corresponding 
differential cross section diverges as either axis of the plane is approached and at the origin, i.e. when the emission 
becomes soft or collinear. 

The PSS methods regulate the singularities by partitioning the phase space into a region of resolved emissions, and 
a region of unresolved soft and collinear emissions. The resolved part is integrated numerically. The contribution from 
the unresolved soft and collinear emissions is calculated analytically and included with the n-body squared matrix 
element such that the net result is finite, though not necessarily positive. For the case of Smin-slicing, the boundary 
of the unresolved region is defined by a single parameter with dimension energy squared. An emission is considered 
to be unresolved anytime the invariant mass squared of any parton pair is less than the Smin resolution parameter 

\sij\ < Sinin (uurcsolvcd rcgion) (2) 



where the partons i, j may be either initial or final state. 

The cross section for a particular n-body configuration $„, integrated over the entire ii vs. i plane of Figure ^ is a 
constant, schematically 

'T"($„,S,ni„) + / (T"+l($„,$ + l)d$+l = Const($„) (3) 

and is independent of the Smin choice. For a complete description of Smin-slicing refer to [ p7| , p^ . 

There exists a boundary in the unresolved region of t vs. u space, referred to here as Szcro, inside which the sum 
of the n-body and (n-l-l)-body contributions is zero. Knowing the location of the boundary, one could calculate the 
NLO cross section and distributions by sampling only the (n-l-l)-body phase space, restricted to that region which 
lies above the boundary (resolved partons). Thus cr"(<l>n, Szoro) = and the constant of Eq. ^ is 

f a"+i($„,$+i)rf$+i = Const($„). (4) 

The analytic expression for cr"($n,Smin) is given in Ref. JTH , and the analytic expression for Szcro is derived in the 
appendix of the present paper. 

A variant of the idea was originally proposed by Baer and Reno in Ref. p9[ | , who approximated the Szcro boundary 
as constant and evaluated it by trial and error for single vector-boson production in hadronic collisions using the two 
parameter PSS method. However, the location of the boundary varies event-by-event with the n-body kinematics. 
This was demonstrated by Potter [p| , who formulated techniques for evaluating the Szero boundary dynamically event- 
by-event. The idea has been implemented for jet production in deep inelastic scattering in Ref. pO| , and good 
agreement is found with the Hera data. 



III. THE PHASE SPACE VETO METHOD 

To illustrate the phase space veto ($-space Veto) method for hadronic collisions, the process pp^ ^ Z'^ + X ^ 
-\-X is chosen. As for the PSS methods, the phase space is divided into two distinct volumes. The n-body volume 
encompasses the phase space with no resolved emission (pp' — s- Z*^ — > I'^l^ kinematics), while the (n-l-l)-body volume 
describes the phase space with an extra parton in the final state, pp^ — > Z'^j l^l^j, where j denotes a gluon or 
(anti)quark. For the PSS method, the two volumes would be integrated separately using numerical techniques, and 
then added together. 

For the <i>-space Veto method, the integration is organized differently. Only the (n+l)-body volume is integrated, 
and the n-body matrix elements are used to test on which side of the Szoro boundary each phase space point lies. An 
event candidate sampled in the (unrestricted) (n-l-l)-body phase space represents a point in the u vs. t plane shown 
in Figure |^. If the point lies below the Szoro boundary, the event is veto-ed. If it lies above the boundary, it is assigned 
the event weight from the (n-l-l)-body differential cross section. Since the location of the Szcro boundary depends on 
both the factorization and renormalization scales, the reduced scale dependence of the NLO calculation is maintained. 

There are at least two possibilities for determining on which side of the boundary a phase space point lies: 

1 . the location of the Szoro boundary can be calculated analytically. The n-body cross section is a quadratic equation 
in InSniiniWith the smaller of the two roots corresponding to the correct solution. This is the method proposed 
in Refs. [py20[| , where the Szom equations for single jet production in electron-proton scattering are derived. In 
the appendix of the present paper, the corresponding Szcro equations for pp ' —> at NLO are derived. 

2. without knowing the location of the Szoro boundary in the unresolved region explicitly, it is possible to test 
on which side of the boundary a phase space point lies by projecting the (n-f l)-body kinematics onto n-body 
kinematics and simply evaluating the sign of the n-body matrix element with the Smin boundary adjusted to 
sit on top of the point in the u vs. t plane. It is not necessary to keep track of Jacobians from the projection 
nor overall normalization factors, since only the sign of the matrix element is of interest. One must be careful 
because well above the Szoro boundary (and after the Smin-slicing approximation has broken down) the n-body 
cross section turns negative once again (corresponding to the second solution of the quadratic equation in InSmin, 
discussed above). In practice this happens only at large (« lO^GeV^) values of Smin ^ lOOGeV^ (see Figure ||). 



This strategy is simple to implement, and works for Smin-slicing, one parameter PSS, and two parameter PSS.|^ In 
this manner, processes which have already been coded as a numerical integration using one of the PSS methods 
can be re-cast as event generators with minimal effort. 

Regardless of which of the above techniques is chosen, it is necessary to project 7 dimension (n-l-l)-body kinematics 
onto the 4 dimension n-body ones. This is accomplished by requiring the lepton-pair mass Mj+i- , and rapidity Yj+j- , 
to remain unchanged in the projection. To perform the projection, the center of mass frame lepton momenta are 
boosted into the vector-boson rest frame (which is the 'new center of mass frame'), and then boosted longitudinally 
such that the pair regains their original rapidity, Yi+i- . 

In Figure |4| the Szcro boundary for the Tevatron and LHC collider energies are shown as a function of the lepton- 
pair rapidity for several parton center of mass choices. The dependence of the Szero boundary on the choice of 
renormalization and factorization scales is shown in Figure 

Though the Szcro boundary always exists, there is no guarantee that the boundary lies within the region of validity 
for the PSS methods. This has not been a problem for the limited set of processes to which the method has been 
applied. However, a hybrid of the PSS and subtraction methods has been proposed in Ref. ||^ to deal with the 
situation, should the need arise. 

For each phase space sample in the above algorithm, both the (n+l)-body and n-body matrix elements are evaluated. 
This means that the event generation will be slower than that for tree level events by the amount of computer time 
it takes to evaluate the n-body matrix elements which are used to perform the veto. Though this appears to be the 
minimal computation necessary for performing a calculation which incorporates the full NLO information, this is not 
the case. There are ways in which the performance, in terms of computational time, can be improved: 

• upper and lower limits on Szcro can be evaluated (see Figure For phase space points which lie outside of these 
limits, the n-body matrix element need not be evaluated to determine whether or not the point is veto-ed. For 
Tevatron energy, Szcro ranges from about 1 GeV^ to about 100 GeV^. 

• since event generation is normally implemented using the hit-and-miss (i.e. acceptance/rejection) Monte Carlo 
technique, the majority of event candidates will be rejected anyway. The $-space Veto need be applied only to 
those event candidates which are accepted (or whenever an event candidate violates the maximum event weight 
against which the acceptance/rejection is taking place). Since the efficiency of event generators is typically 
about 25% or lower, it means that the n-body matrix element needs to be evaluated rarely. Further, when the 
event candidates are sampled from an adaptive integration grid (such as for the implementation presented in 
this paper) , the adaptive integration will "learn" the location of the boundary, and will bias the sampling away 
from the region below the boundary. 



IV. NUMERICAL RESULTS 

The pp' ^ Z'^ + X ^ + X event generator is implemented using the squared matrix elements of Ref. ||2^ cast 
into the Smin-slicing method [ psf (which employs special 'crossed' structure functions). The matrix elements include 
both the and 7* diagrams with decay to massless leptons, such that the branching ratio to one lepton flavor is 
automatically included. This means finite width effects, lepton decay correlations, and forward-backward asymmetries 
are everywhere taken into account. The generator is written in C-I--I- using modern object-oriented design patterns. A 
new prototype C-I-+ version of the Bases/Spring program is used for adaptive integration and event generation. 
Special care has been taken to make the program user friendly, and it is available upon request from the author. 

All of the distributions and cross sections presented in this paper are for pp collisions at 2 TeV (Tevatron Run II) or 
pp collisions at 14 TeV (LHC), with the lepton-pair mass restricted to the range 66-116 GeV and decaying to e~e^ . 
CTEQ3M parton density functions are used (chosen because the 'crossed' versions of the structure functions p^ ] 
are readily available, though in principle they can be tabulated for any structure function). For all calculations the 
renormalization and factorization scales have been set equal to the vector-boson mass, = fJ-p = Afi+i- , and the MS 



^ the two parameter PSS method must first be expressed in terms of a single parameter, for example by defining 5c ~ O.ISb 
where 5c and 5b are the coUinear and soft parameters of the method. The author has tested this for pp ^ Z'^ + X at NLO 
and found good agreement both with the unaltered two parameter PSS method, and also with the "l>-space Veto distributions 
presented in this paper. 

^ Hence the vector-boson is denoted by Z" , even though the 7* contribution is included. 



factorization scheme is used. The input parameters are chosen to coincide with those in PYTHIA 6.200: the Z° mass 
and width are Mz — 91.188 GeV and Tz = 2.47813 GeV, the electroweak mixing angle is sin^ 9w = 0.232, and the 
electroweak couphng is aQED(Afz) = 1/128.8. The two-loop MS expression for as is used with A^''^^ = 0.239 GeV. 
Using these input parameters, the <i>-space Veto event generator predicts 239.7 ±0.6 pb for the inclusive cross section 
at Tevatron Run II. 

In Figure ^ the inclusive cross section prediction from the <I>-space Veto event generator is compared to the pre- 
dictions from the Smin-slicing calculation using several choices of the s,„i,i parameter. The results are consistent, 
indicating the Szoro boundary lies within the region where the Smin-slicing approximation is valid. 

In Figure distributions produced with the ^-space Veto event generator are compared to those derived from nu- 
merical integrations using Smin-slicing. The $-space Veto method faithfully reproduces the NLO transverse momentum 
of the electron. The transverse momentum of the vector-boson also agrees well with the Smin-slicing everywhere that 
the NLO calculation is valid. 

In the small region, multiple gluon emission becomes important and fixed order perturbation theory is unreliable. 
This is evident in the inset of Figure ^ In this region the results depend on the specific choice of the Smin parameter. 
This is also the region where the <I>-space Veto method becomes unreliable because the minimum jet scale is coupled 
to the n-body kinematics. This effect is visible in Figure ||, where the kinematics of the veto-ed event candidates from 
the <I>-space Veto method for a typical event generation run are plotted in the vs. Smin plane. The largest P^ 
of a veto-ed candidate event is 5.5 GeV, indicating the NLO calculation is unable to provide a useful prediction in 
the region below « 5.5 GeV. It is worth stressing that this does not make the $-space Veto method less useful than 
Smin-slicing since any NLO calculation is unreliable here. This is the region where the distributions are better modeled 
with the parton shower, and a suitable treatment which removes this minimum jet scale coupling will be provided 
in the next section. The Szoro boundary represents a lower limit to the usefulness of our fixed order perturbative 
approximation. As such, Szcro is a useful concept as a qualitative measurement of the frontier of the validity of our 
perturbative calculation. 

In Figure |^ the factorization and renormalization scale dependence of the transverse momentum of the electron 
and vector-boson distributions are shown using the <I>-space Veto method for Tevatron and LHC energies. The scale 
dependence is identical to that from the Smin-slicing method, because the Szero boundary encodes information about 
the scale choices (Figure The change in the distributions resulting from the variation of the scales is an indication 
of the theoretical error arising from neglected higher order terms. The importance of the reduced scale dependence is 
demonstrated in Figure |l^, where the variation in the prediction at Born level and at NLO of the transverse momentum 
of the electron for pp ^ Z'~' + X ^ e+e~ + X a.t LHC energy is shown. The comparison is restricted to that region 
where the Born level prediction is meaningful. The same comparison is shown for the lepton-pair mass distribution in 
Figure |ll|. The change in the scale at Born level results in more than a 25% variation in the distributions, whereas the 
prediction from the NLO 'I'-space Veto generator reduces this variation to about 7%. The scale dependence arising 
in predictions from event generators which use leading order subprocesses (like PYTHIA, HERWIG, and ISAJET) will 
resemble that of the Born level prediction. 



V. SHOWER EVOLUTION 



At the present stage, each event consists of the vector-boson decay products and exactly one colored emission in the 
final state. The energy scale of the emission is at least y'Szero. Unweighted events are provided by the Bases/Spring 
algorithm, and the normalization is NLO. A coupling between the minimum emission scale y'Szero and the kinematic 
configuration exists in the very small P^ region. 

The next step is a consistent interface to a parton shower algorithm. The goal is to have the parton shower dominate 
the prediction in the soft/coUinear region (in particular, it should preserve the parton shower's prediction of Sudakov 
suppression ||25[|), and the first order tree level diagrams dominate in the region of hard well separated partons. This 
does not compromise the integrity of the prediction, it merely highlights that different approaches are well-suited to 
different regions. 

To accomplish this, a parameter Sp.s. is introduced to partition a region of t vs. u space which is exclusively the 
domain of the parton shower. This parameter may be thought of as separating the fixed order regime from the 
all-orders parton shower region, in the same way that a 0(1 GeV) parameter in the showering and hadronization 
programs defines the scale at which the parton shower is terminated, and the simulation turns to the non-perturbative 
hadronization model for a description of the physics. This partition is shown in Figure |l^. Events which lie below 
the Sp.s. boundary are first projected onto n-body kinematics (i.e. the point in t vs. u space is moved to the origin) 
and the parton shower is allowed to evolve the event out into the plane. The projection is performed keeping the 
lepton-pair mass and rapidity fixed, exactly as described in Sec. III. The parton shower is invoked with the scale set 



to ^Sp.s., which ensures the evolution does not move the event out into a region of phase space which has aheady 
been counted using the first order tree level matrix elements. 

A reasonable choice for the ^sp.s. parameter is a few times the minimum jet scale, ^^Szero- This ensures the first 
order tree level matrix element is reliable above the sp.s. boundary. The distributions have very little sensitivity to 
the choice of sp.s.. 

For events which lie above the sp.s. region, the parton shower is also invoked, this time with a scale equal to the 
minimum invariant mass of any parton-pair 

parton shower scale — minimum [ Q^g, Qqg-, Qqg] (5) 

which ensures no double counting can occur. 

The effect of the projection and subsequent parton showering is shown in Figure ^ Initially the distributions are 
provided by the <i>-space Veto, solid line. The projection is applied to events which sit below the sp.s. boundary, which 
effects only the small region, and is shown as a dashed line and does not correspond to anything physical. Finally 
the parton shower is applied (dotted line), and has the largest effect on those events which have been projected. 

In Figure |l|the $ -space Veto event distributions (including parton shower evolution) are shown for several choices 
of the Sp.s. parameter. The dependence on the Sp.g. parameter choice is small, indicating discontinuities which might 
exist at the sp.s. boundary are also small. 

For the distributions presented here, events from the <I>-space Veto generator have been evolved with the PYTHIA 



6.200 parton shower. PYTHIA is attached using the HepUP interface |26|, which is a generic standard for the com- 
munication between event generators. Having evolved the events through the parton shower, PYTHIA provides other 
features of the event structure such as hadronization, resonance decays, beam remnants, and multiple interactions. 
The showered event distributions presented in this paper include all of these features. The use of the HepUP interface 
allows for the parton shower program to be easily interchanged. The choice of PYTHIA is arbitrary, there is nothing 
which precludes the use of any other showering program.^ 

The full event generator is now complete. Adaptive integration and phase space generation is provided by 
BasesSpring. The event weights are evaluated using the $-space Veto method, which discards those event candi- 
dates lying below the Szoro boundary. When the program is executed, the phase space is first mapped onto a grid 
using an initialization pass with the adaptive integration (performed by the 'Bases' part of the BasesSpring package). 
The 'Spring' part of the BasesSpring package then provides unweighted events, by sampling candidate events from 
the adaptive integration grids and accepting events according to the differential cross section using the acceptance- 
rejection algorithm. After removing the emission from those events which are soft or collinear (as defined by the sp.g. 
parameter), the events are transferred to the PYTHIA package using the HepUP interface. PYTHIA performs the parton 
shower, and subsequent event evolution including hadronization, etc. 

While the end result in terms of physics does not differ significantly from that obtained by the author in Ref. j|] 
for W Z production, the method presented here is simpler, easier to implement, faster in terms of computer time, and 
may be generalized to a broad range of processes. Improved methods for invoking the parton shower from parton 
level event configurations are being developed [^,^, and are suitable for application to the $-space Veto events. ^ 

A comparison of the computer time for generating the events is presented in Table |. The processing time per 
event and generation efficiency (percentage of candidate weighted events which are accepted in the event generation 
algorithm) for PYTHIA and the $-space Veto are similar, indicating the ^-space Veto method is successful in encoding 
the extra NLO information without affecting the overall time-performance of event generation. 



In Figure 15 the $-space Veto distributions (solid line, includes evolution with the PYTHIA showering and hadroniza- 
tion package) are compared to the predictions from PYTHIA. In PYTHIA there are two strategies implemented for single 
vector-boson production. For both strategies the hard subprocess is chosen according to the Born level matrix ele- 
ment, such that the normalization is always leading order. For the "old" PYTHIA implementation of the process, the 
event is then evolved with the standard parton shower beginning at a scale equal to the vector-boson mass. For the 
new "matrix element (M.E.) corrected" implementation of the process, the shower is initiated at a scale equal to 
the machine energy and is corrected according to the Z^-fjet first order tree level matrix element, which results in a 



^ For the case of the HERWIG parton shower, there is a region or "dead zone" in the i vs. u plane of Fig. g where emissions 
never occur. The boundary of the dead zone is a natural choice for the partition which separates the parton shower region 
from the region populated directly by the first order matrix element when using HERWIG. This is the prescription employed in 
Ref. for "hard matrix element corrections" to single vector-boson production. 

* The apacic++ showering program employed in Ref. does not yet include initial state showers, but an implementation 
is expected soon. 



considerable improvement of the high region modehng. The virtual one- loop contribution is not included anywhere 
in the PYTHIA implementations. The dotted line in Figure |l^ is from the standard PYTHIA process, and the dashed line 
is from the M.E. PYTHIA process. The ^-space Veto distribution and M.E. corrected PYTHIA shapes are rather similar, 
indicating the matrix element corrections in PYTHIA are having the desired effect. The $-space Veto distributions 
have the advantage of NLO normalization and a reduced dependence on the factorization and renormalization scales. 



VI. CONCLUSIONS 



The $-space Veto method for organizing NLO calculations into event generators is demonstrated for production 
in hadronic collisions. The method is based loosely on the ideas proposed by Potter for deep inelastic scattering [||. 
The primary motivation for the method is to move numerical NLO calculations beyond the status of "event integrators" 
to "event generators" , making them suitable for interface to showering and hadronization programs and subsequent 
detector simulation. 

The general features of the $-space Veto method are: 

• event weights are positive definite, meaning the standard methods for event generation can be applied, providing 
a prediction which is well suited for experimental applications. 

• in the soft/coUinear region, the results are dominated by the parton shower. In particular the low region 
exhibits Sudakov suppression. 

• in the region of hard well separated partons, the distributions are dominated by the first order matrix element. 

• the normalization is NLO and the reduced scale dependence afforded by the NLO calculation is maintained. 

The method has been implemented as an event generator (available from the author) for pp^ — > /^* + X — > 
I'^l^ + X, with showering and hadronization provided by the PYTHIA package. 
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APPENDIX: THE s.oro FUNCTION FOR pp ^ + X AT NLO 

The differential cross section for the n-body contribution to pp' Z^ + X evaluated in the MS scheme, integrated 
over unresolved emissions out to a scale of Smin, and neglecting overall factors, is a quadratic equation in InSmin 

(6) 
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where the sum runs over all flavors of initial state (anti)quarks, Nc = 3 is the number of quark colors, y/§ is 
the vector-boson mass, A^Born ^^e Born level matrix element for q'^q'^ — > Z^, f p_i--) {x, ^J■F) is the parton density 



function evaluated at Bjorken momentum fraction x and factorization scale /if, the renormalization scale is (often 
A'F = M-R = is chosen), and 



P-i q ■ 



are the crossed structure functions presented in Eq. 3.37 of Ref. jlj]. 

The solution for Szoio corresponds to the smaller of the two roots of Eq. ^ 



with 



minimum 



s X exp 



-b - vb2^^4ac 
2a 



, s X exp 



-b + Vb2 - 4ac 
2a 



b = E,, ~^■m^\^\MtJ'.fp^^^ix,,^,F)fp.^,i^2,^^F) 

+ ^#^A^c|A^Lnl' [fp^^^ixu^^F)Ap,^,{x2,^iF) + Ap,^,{xu fip) fp^^j {X2, f^p)] 



Eij 



asiPR) Nl-l f^TT^ . 1 _ _9 
2tt Nc 



+ ^^Nc\Mtj' [fp.-.^{xl,fiF){Ap,^,{x2,HF)ln^ + (x^, mf)) 



(7) 



(8) 



(9) 



+ [Ap,^,{xi,^iF)lll-4- + Bf°_^,^{xi,^lF)\ fp^^j{x2,^iF) 
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real emission 
FIG. 1. Feynman graphs contributing to pp 
vector-boson decay products are not shown. 



+ X at NLO. The wavy line represents either a or 7*, and the 




FIG. 2. A projection of the pp Z^j phase space onto the u vs. i plane is shown, where u = (p2 — pj)^ = —Q2j 
t — {pi — PjY — ~Qij^ P2, Pj are the momenta of the forward colliding parton, backward colliding parton, and real 

emission. The area above (below) the Smin boundary is the region of resolved (unresolved) real emissions. When Smin = Szcro, 
it denotes the boundary defining the region inside of which the the n-body and (n+l)-body contributions sum to zero (i.e. the 
cross section integrated over the unresolved region is zero). 
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FIG. 3. The two roots of the quadratic n-body differential cross section presented in Eq. ^ are plotted as a function of the 
lepton-pair rapidity, evaluated at parton center-of-mass energy equal to the mass for pp collisions at 2 TeV (Tevatron, top) 
and for pp collisions at 14 TeV (LHC, bottom). The smaller solution is the Szcro function of interest, the larger solution should 
not be interpreted physically. 
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FIG. 4. The dependence of Szcro as a function of lepton-pair rapidity at several choices of parton the center-of-mass energy 
Q is shown for the pp colhsions at 2 TeV (Tevatron, top) and for pp coUisions at 14 TeV (LHC, bottom). The Szero function 
does not depend strongly on the vector-boson decay angles. 
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FIG. 5. The scale variation of the Szcro function evaluated at parton center-of-mass energy equal to the mass is shown for 
pp collisions at 2 TeV (Tevatron, top) and for pp collisions at 14 TeV (LHC, bottom). The Szcro function encodes information 
about the factorization and renormalization scale choices into the ^-space Veto method, preserving the NLO calculation's 
reduced scale dependence. 
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FIG. 6. The inclusive NLO cross section for pp X ^ e^e^ + X reactions at 2 TeV with the lepton-pair mass 

restricted to 66-116 GeV is shown as a function of the Smin parameter for the Smin-sHcing method. The cross section calculated 
using the *l>-space Veto event generator is superimposed and is in good agreement. The Born- level cross section is also shown. 
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FIG. 7. The transverse mornontum of the electron (top) and vector-boson (bottom) are shown for the process 
pp — > Z° + — » e^e~ + X at 2 TeV with the lepton-pair mass restricted to 66-116 GeV (no parton showering is used). 
Distributions derived from numerical integrations at NLO using Smin-slicing for various choices of the Smin parameter are com- 
pared to the distributions from the NLO $-space Veto event generator. Agreement is excellent everywhere, except in the low 
Pz region (inset) where fixed order perturbative QCD is unreliable. The Born level prediction is also super-imposed for the 
Pj- distribution (top). The Born level prediction for the vector-boson transverse momentum is a delta function at P[ = 0. 
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FIG. 8. The kinematics of $-space Veto method event candidates which have been veto-ed because they lie below the Sz 



boundary are plotted in the Pz vs. 



plane. The largest of a veto-ed candidate event is 5.5 GeV. The process is 



pp — > Z -h X — » -\- X aX 2 TeV with the lepton-pair mass restricted to 66-116 GeV. 
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FIG. 9. The transverse momentum of the electron (top) and vector-boson (bottom) are shown for the process 
pp — > Z° + X — > e"'"e~ + X at NLO using the $-space Veto method (no parton showering is used) for different choices of 
the renormalization and factorization scales. The spread in the distributions is an indication of the theoretical error from 
neglected higher order terms. The distributions on the left arc for 2 TcV pp collisions at the Tevatron, and the ones on the 
right are for 14 TeV pp collisions at the LHC. The lepton-pair mass is restricted to 66-116 GeV, and the three curves use the 
same event sample. 
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FIG. 10. The reduced scale dependence of the NLO 3>-space Veto calculation as compared to the Born one is demonstrated. 

The transverse momentum of the electron (top) for the process pj) Z'^ + X ^ e'^e~ -\- X with the lepton-pair mass restricted 
to 66-116 GeV at LHC energy (14 TeV) is shown (the vector-boson transverse momentum is not shown because the Born level 
calculation does not provide a prediction for P^). The renormalization and factorization scales are varied by a factor two in 
the Born level calculation (left) and the NLO $-space Veto calculation (right). The percent variation of the distributions is 
shown at bottom. The variation is about a factor 4 smaller for the NLO $-space Veto. The effect is smaller at Tevatron energy 
in this region. 
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FIG. 11. The reduced scale dependence of the NLO $-space Veto calculation as compared to the Born one is demonstrated. 
The lepton-pair mass in the vicinity of the Z'^ resonance is shown (top) for the process pp Z'^ + X ^ e^e~ + X at LHC 
energy (14 TeV). The renormalization and factorization scales are varied by a factor two in the Born level calculation (left) 
and the NLO 3>-space Veto calculation (right). The percent variation of the distributions is shown at bottom. The variation is 
about a factor 3 smaller for the NLO $-space Veto. The effect is smaller at Tevatron energy. 




FIG. 12. A projection of the pp Z'^j phase space onto the u vs. i plane is shown, where ii — (p2 — Pj)'^ = — Qij and 
t = (pi — Pj)^ = ~Qij^ a'^d Pi' P2' Pj ^■'s ^lis momenta of the forward coUiding parton, backward coUiding parton, and the 
hardest emission. Events in the region of hard weU separated partons are sampled with the first order matrix element, then 
evolved further by the paxton shower. Events in the region between the Szero and sp.s. boundaries are projected onto n-body 
kinematics (i.e. onto the origin of the plane) and then evolved with the parton shower to a point which may lie anywhere below 
the Sp.s. boundary. The region below Szero is never sampled, but may be reached by the projected and showered events. 
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FIG. 13. The Pj distribution is shown after different stages of the event generation for pp ^ Z" + X ^ e^e~ +X at 2 TeV 
with the Icpton-pair mass restricted to 66-116 GcV. The solid line is the <I>-spacc Veto NLO distribution without any projection 
or parton shower. The dashed line is the (nonphysical) distribution for the same event sample, after applying the projection 
with ^sp.s. = 25 GeV. The dotted line is the distribution after subsequent evolution through the showering and hadronization 
program. 
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FIG. 14. The effect of the parton shower on the $-space Veto distributions is shown for several choices of the sp.s. parameter 
which partitions the phase space into the region populated by the paxton shower, and the region populated directly by the first 
order matrix elements. The transverse momentum of the electron (top) and vector-boson (bottom) are plotted for the process 
pp ^ Z'^ + X ^ e^e~ + X at 2 TcV with the lepton-pair mass restricted to 66-116 GeV. There is very little dependence on 
the specific choice of the sp.s. parameter. 
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FIG. 15. Distributions for the process pp —* Z'^ -\- X —> e^e~ + X at 2 TeV from the $-space Veto event generator (sohd 
hne, includes evolution through the PYTHIA shower and hadronization program) are compared with the PYTHIA internal process 
distributions. The dashed line is the "matrix element corrected" PYTHIA prediction and the dotted line is the "old" (no matrix 
element corrections) PYTHIA prediction. The lepton-pair mass is restricted to 66-116 GeV. 

TABLE I. A comparison of computer processing time for the $-space Veto method and for PYTHIA. In each case 10000 
events are generated for the process pp ^ Z'^ -\- X ^ e^e~ + X eX 2 TeV and the events are evolved through the PYTHIA 
shower and hadronization program. The $-space Veto event generator requires some time to initiaUze the Bases/Spring grids 
(i.e. "learn the phase space"), whereas no initialization time is required for PYTHIA processes. The processing time per event 
and efficiency are similar. The computer is a 650 MHz Pentium III. 



Method 


Time for Grid Initialization 


Time for 10000 Events 


Efficiency 


$-space Veto 


14.0 s 


70.3 s 


28% 


PYTHIA 




68.6 s 


27% 



